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ABSTRACT: Recently Callan, Giddings, Harvey and the author derived a set of one- 
loop semiclassical equations describing black hole formation/evaporation in two-dimensional 
dilaton gravity conformally coupled to iV scalar fields. These equations were subsequently 
used to show that an incoming matter wave develops a black hole type singularity at a 
critical value <p cr of the dilaton field. In this paper a modification to these equations arising 
from the Fadeev-Popov determinant is considered and shown to have dramatic effects for 
iV < 24, in which case (j) cr becomes complex. The N < 24 equations are solved along 
the leading edge of an incoming matter shock wave and found to be non-singular. The 
shock wave arrives at future null infinity in a zero energy state, gravitationally cloaked by 
negative energy Hawking radiation. Static black hole solutions supported by a radiation 
bath are also studied. The interior of the event horizon is found to be non-singular and 
asymptotic to deSitter space for iV < 24, at least for sufficiently small mass. It is noted 
that the one-loop approximation is not justified by a small parameter for small N. However 
an alternate theory (with different matter content) is found for which the same equations 
arise to leading order in an adjustable small parameter. 



In recent work [1] it was proposed that two-dimensional dilaton gravity coupled to con- 
formal matter is a useful and simple model in which progress might be made in unravelling 
the mysteries associated with black hole evaporation [2]. It was shown that the process 
of black hole formation/evaporation, in an approximation which includes one- loop matter 
effects and treats gravity semiclassically, can be described by a set of partial differential 
equations which incorporate the backreaction of Hawking radiation on the geometry. It 
was further pointed out that this approximation is formally exact in a limit as the number 
N of matter fields is taken to infinity. 

It was subsequently shown [3,4] that in the large N approximation a collapsing matter 
wave still forms a black hole containing a singularity*. This singularity no longer occurs at 
the value <fi = oo of the dilaton as in the classical theory, but rather moves up to the finite 
value (j) cr = —^injrj. The black hole then evaporates, presumably leaving a massless, 
singular "remnant" [5,6,7]. 

In this paper we will consider the equations describing dilaton gravity coupled to N 
conformal scalars in the one-loop approximation for finite N. These equations differ from 
those derived in [1] by the addition of terms arising from the gravity-ghost measure which 
are negligible for N^oo. We shall see that for N < 24, these terms remove the singularity 
found in [3,4]. 

Some evidence consistent with the absence of other types of singularities is presented, 
but the equations are sufficiently complex that the question is not settled here. We hope 
to analyze the problem numerically in the near future [8]. 

We begin with a discussion of the gauge fixing and quantization of pure dilaton gravity 



* Or more precisely, the fields become so large that the large- A?" approximation breaks down. 
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(diet) 



where g and (p are the metric and dilaton fields respectively and A 2 is a cosmological 
constant.* This is a theory with no local degrees of freedom: the 3 + 1 fields in g and 
may be eliminated by two gauge conditions and two constraints. There are however a one 
parameter family of classical black hole solutions labelled by the black hole mass [10]. We 
wish to gauge fix (diet) to conformal gauge 

9++ = 9 — = 0, (gfer) 

where a ± = r±a. In so doing the action will be shifted by the usual logarithm of the 
Fadeev-Popov ghost determinant. This term may be expressed in a covariant notation as 
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J d 2 a^gErT l R. (fpet) 



However an ambiguity in this procedure arises in the present context. There is a family of 
metrics g a given by 

9a = e^g (glph) 

any of which might be used to construct Spp in (fpet). Although we have chosen go in 
order to write down (diet), it is not especially preferred. Indeed the difference between 
two choices of metric is given by 

S F p(9a) - S FP (g) = J^J d 2 a^g~(R + oo0)0. (delfp) 

This is a local expression which might have been added to the action (diet) either in the 
first place or as a finite counterterm during one-loop renormalization. Thus there is no 
right or wrong choice of metric in (fpet): different choices simply correspond to different 
theories. One must choose a theory which contains the physical phenomena one wishes to 
investigate. 



* A related discussion of this quantization with similar conclusions has been given in [9] . 
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In fact g is not a good choice of metric to use for denning the ghost measures-it leads 
to a sick theory. Using the metric g in (fpct) means that the Fadeev-Popov b — c ghosts 
couple to the geometry in the same way as the conformal / fields described in [1]. It 
immediately follows that black holes will grow in mass by Hawking radiation of negative 
energy ghosts.* This is clearly nonsense. 

This problem is avoided by defining the ghost measure with the alternate metric 

9 = e" 2( V (ghat) 

This metric turns out to be flat for all classical solutions of (diet). Black holes will therefore 
not radiate ghosts to leading order. 

As is familiar in Liouville gravity, there is an additional term of the form (fpct) arising 
from the dependence of the p, measures on the metric. Again there is an ambiguity in 
these measures. Since p, are not local, propagating degrees of freedom, it is natural to 
demand that there be no Hawking radiation in these modes. This is accomplished by using 
the metric g to define their measures as well, which changes the 13 in equation (fpct) to a 
12. This definition ensures that there is a stable black hole solution of the quantum theory 
for each value of the mass M. 

The gauge fixed action, including all the measure terms, is then 

Sq+S m = ^J d 2 a[{e- 2 ^{2d + d-p - 49+09-0+ X 2 e 2p ) 

+ 2d+(p-0)cL(p-0))], (rhct) 

* Although we cannot build detectors to see the ghosts directly, we can still observe their 
effects on the geometry. 
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while the stress tensor is** 

Tl + +Tf + = e- 2<t) {4d + (i)d + p-2dX<f)) + 2(d+p - <9+0) 2 - 2d\p+2d\<j) + t+. (rht) 

where, as explained in [1], t+(a + ) is determined by boundary conditions. 

In order to study the problem of black hole evaporation we must complicate the theory 
by adding matter fields with local degrees of freedom. Following [1] we add N confor- 
mally coupled scalar matter fields /j. Again an ambiguity in defining the / measure arises. 
However this time we do not wish to use the metric g. This leads to a presumably sen- 
sible theory which does not contain the phenomena we wish to study: black holes do not 
Hawking radiate. Indeed this theory in a sense does not even contain black holes, since 
the matter sees only the flat metric g. As explained in [1], if we instead use g to define 
the / measure Hawking radiation of / particles indeed occurs, and closely resembles the 
four-dimensional phenomena. One thereby arrives at the final action 

S = IJ d 2 a[e- 2(t> (2d+d-p - 4<9+0d-</> + \ 2 e 2p ) 

N 1 N 

- - d+pd-p + - J2d+fid-fi+2d + (p - 0)d_(P " 0)]- (ghct) 

i=l 

The constraints will be discussed shortly. 

The quantum theory is described by functional integration with the "naive" measure 
weighted by S. In the large- N limit, all terms are of order iV (after shifting 0) except the 
last one, which is order one and may therefore be dropped. In [1] it was argued that these 
order N terms may be treated as a quantum effective action which describes the process 
of black hole formation and evaporation, with the modifications of the gravitational action 

accounting for the stress-energy carried by Hawking radiation. 

** Ignoring the cosmo logical constant term, (rhct) becomes a free theory in terms of the 

variables v = e~ 2< ^ and w = p — 4>. It is then easy to check that the stress tensor has 
c = 26, as required by coordinate invariance. It also easily follows that the cosmological 
constant operator e w is dimension (1,1) (with no renormalization of the exponent). Thus 
there are many similarities with Liouville theory, and the methods developed there may 
be useful in the present context. 
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However, the large- A approximation is not necessarily a reliable barometer of finite- A 
physics, particularly for A < 24. One way to see this is from the behavior of the p— 4> kinetic 
operator /C at large positive and negative <fi. At large negative 0, the theory is essentially 
classical and /C has one positive and one negative eigenvalue. At large positive and finite 
A, the classical action may be treated as a perturbation about the free measure-induced 
terms, and /C still has one positive and one negative eigenvalue*. In the large-A limit 
however, there are two negative eigenvalues for large 4> and consequently a zero eigenvalue 
at an intermediate value of 4>. This zero eigenvalue leads to singular behavior in the large- 
A limit [3,4] which may not be present for small A. Thus other methods should be found 
for analyzing the theory at small A. 

In this paper the action (ghct) including the last term (and corresponding modification 
of the constraints) will be treated as a quantum effective action for finite A. This amounts 
to a one-loop semiclassical approximation. For small A, there is no obvious small parame- 
ter which justifies this approximation. The loop expansion may break down when e 2( ^ gets 
large, and we cannot be confident that our conclusions are qualitatively correct. Neverthe- 
less, we shall find in the one-loop approximation that the behavior of the theory changes 
dramatically at A = 24, and we hope that the one-loop semiclassical approximation is at 
least qualitatively correct in the A < 24 regime. 

While treating (ghct) semiclassically has not been justified as a systematic approxima- 
tion to dilaton gravity coupled to A scalar fields for small A, it can be formally justified 
as a systematic approximation to dilaton gravity coupled to a different matter system: let 
there be AM scalar / fields, and include an additional c = — 24M conformal matter sector 
with measure defined with respect to g. After a shift of 0, one recovers an action of the 
form (ghct) multiplied by M. One then expects a semiclassical treatment to be valid for 

* Although for A > 24 there is a region whose size grows with A in which there are two 
negative eigenvalues. 
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large M. The following analysis may be taken to apply to this system. 

We now proceed to analyze the dynamics following from (ghct) . The p and equations 
may be cast in the useful form 

8Pd+d-(p = -P'(4d+(f)d-(j) + X 2 e 2p ), (pheq) 
2Pd+d-p = e~ A(t, {Ad + (j)d-<P + X 2 e 2p ), (rheq) 

where 

12 + 24' 

The ++ constraint equation is 

1 N 

T++ = e-^^d+^d+p - 2^0)+- J^d+fid+fi 

i=l 

-^(d+pd+p - d 2 + p)+2(d+(p - <P)d+(p - 0) - d\ (p - 0)) + t+ = 0(+ppl) 
and a similar equation holds for T 

The effect of the ghost induced-terms in these equations is immediately evident from 
(pheq). The prefactor P in (pheq) has zeros at 



■m N , / 24 , . 
< r2 *=2A {1 H 1 ~N ) - {ZT0S) 

As pointed out in [3,4], these zeros are very dangerous: because the RHS of (pheq) is 
generically non-zero, d+d-(p is forced to diverge whenever <fi crosses a zero. 

However for iV < 24 there are no real solutions of (zros), and P is a positive definite 
quantity with a minimum at = —^£n ^ = C . f 



f The stability of the zeros (or lack thereof) of P - on which our results strongly depend - 
against higher loop quantum corrections is an important question to which we do not have 
a definitive answer. However it can be easily seen, by considering perturbation theory in 
and e~^, that in the large iV limit (in which the last term in P is neglected), P must 
change sign between weak and strong coupling, and consequently must have at least one 
zero. For iV < 24, it does not change sign, and must therefore have an even number of 
zeros, though we are not sure if that even number is zero in the exact theory. 
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Thus the singularities described in [3,4] do not arise. While we do not know if other 
types of singularities arise in the N < 24 equations, they are clearly far better behaved.* 

These equations can be solved following [3,4] perturbatively about the leading edge of 
an / shock wave incident on the linear dilaton vacuum, as illustrated in Figure 1. The 
linear dilaton vacuum is a solution of (pheq)-(+ppl) given by 

4> = —Act, 

p = 0. (Idvm) 

An / shock wave is defined by 

1 N 
i=l 

Classically, the / stress-energy leads to a black hole of mass M above ctq . The quantum 
equations however exhibit different behavior. Below Oq" one still has (ldvm), and the 
solution above can be computed in a Taylor expansion in (o + — Oq). Defining E(a~) = 
<9+</>(oq~, equation (pheq) becomes a simple equation for E: 

8P8-E = -P'(4d-(pE + \ 2 e 2p ), (sgeq) 

where (ldvm) should be substituted for the values of 0, p along (oq~,ct~). This is easily 
integrated to yield 

E(0 = k i M -A), (sgsl) 

just above the shock wave. The integration constant here is fixed by requiring that asymp- 
totically as cr - ^— oo (on J^) E agree with the classical / shock wave solution. 

For iV < 24, P has no zeros, and E is perfectly finite. (For iV>24, E diverges at a 
finite value of cr~.) As explained in [3], an apparent horizon occurs whenever E vanishes, 



* The weak-coupling singularities of the quantum kink solutions found in [5,7], are presum- 
ably still present. 
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and one may say that an "apparent black hole" has formed. Since the minimum value of 
P (at <]) = 4>c) is 23(1 — 23) there is no apparent horizon for sufficiently weak shock waves, 



i.e., small M. For M = A-y/ ^(1 — ^), E has a double zero where the shock wave crosses 
4> = 4>c, which splits into two apparent horizons (containing a region of trapped pointst) 
as M increases. Since P approaches || as a~^oo (on J+), for M > Ay^ the region of 
trapped points extends all the way from the first apparent horizon up to along the / 
shock wave. 

It is of interest to determine the fate of the apparent black hole-or, equivalently the 
region of trapped points-above the shock wave. This is a difficult problem in general, but 
some progress on it can be made as follows. At the apparent horizon, where d+(f) = 0, 
equation (pheq) reduces to 

P'X 2 e 2 P 

d-(d+<P) = gp- (abc) 

where P' is negative (positive) in the weak (strong) coupling region where </> < C (0 > 0c)- 
It follows that if one moves across an outer (inner) apparent horizon in the direction of 
increasing cr~, one always enters (leaves) the interior of the apparent black hole where 
d+4> > 0. 

Thus if the outer apparent horizon (in the weak coupling region) is followed above 
the shock wave, a + will monotonically increase until (or unless) it meets the line 4> = 4> c 
(which is spacelike inside the apparent horizon). If it does meet this line, the boundary 
(now the inner horizon) must subsequently continue along decreasing a + , back toward the 
shock wave. Thus the apparent black hole ceases to exist for a + greater than the value at 
which the apparent horizon meets 4> c * 



f A trapped point is one point for which <p increases along both outgoing null geodesies. 
This corresponds to the four dimensional definition of a trapped surface when e _2< ^ is 
interpreted as the size of the two spheres [3,4,11]. 

* Though it may (and will for M > Ay^) exist for large values of a~ . 
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It is thus crucial to determine whether or not the apparent horizon actually reaches 
4> c . Following [3], the outer horizon may be parameterized as the curve a~((r + ). Russo, 
Susskind and Thorlacius [3] derive several formulae for a~ from the condition 

d+(j)(a) = d + 2 (f)(a) + ^—d+d-(P(a) = 0. (cdt) 



da+ ^ y ' ^ ' ' da+ 

For the present system these become 

da- __ d\<P 



da + d+d-4> 



X 2 P 
where 

N 



— T® + (aht) 



T ++ = u {d+pd+p ~ d + p) + 2{3+{p ~ 0)<9+(p ~ 0) " 3 + {p ~ 0)) + 1+ {tqu) 

can be thought of as minus the energy in Hawking radiation leaving the apparent black 
hole. The second relation in (aht) implies that if is always negative (i.e. the black hole 
is evaporating) , the outer apparent horizon is timelike and will tend to meet the spacelike 
line 4> = 4> c .] However we have unfortunately been unable to prove that is indeed 
everywhere negative. 

Further progress can be made by considering a very small apparent black hole, for which 
M is just over the threshold for production of an apparent horizon along the shock wave. In 
that case the apparent black hole is formed in a small neighborhood of C , and its evolution 
can be determined by a Taylor expansion about the point (<Tq~, ct^T) = (<7q~, (Jq — j£n ^) 
where the shock wave crosses the line (j) = <fi c - Using the first equation in (aht) and (abc) 
one finds that the trajectory of the horizon is determined by 

da_ w 2(f - Q ^ _ k _ {sbh) 

da + A 3 (cr~ — a c ) a~ — <Jc 



f If goes to zero sufficiently fast the horizon could be asymptotically null and avoid <p c . 
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One then finds for cr+ > <Jq 



cr- = ^2k(a£ - a+) + c 2 + a c (trj) 

where (<7q~, cr^±c) are the initial coordinates of the inner and outer horizons along the 
shock wave. The behavior of the trajectory depends crucially on the sign of d\(f){aQ ,a~) 
or k. According to (trj), if k is positive the apparent black hole shrinks and disappears 
at a + = <Jq + |^7, soon after formation. On the other hand if k is negative the apparent 
black hole will initially grow in size, and perturbation theory about the shock wave can not 
be used to determine its ultimate fate. The sign of k is determined by continuing to one 
higher order the Taylor expansion about <7q used to find £(a~). A tedious computation 
reveals that, rather curiously, k is positive for some values of N and negative for others. 
We do not understand the significance of this. Perhaps the theory depends qualitatively 
on N even within the range < N < 24. 

It is of interest to study the behavior of the quantum stress tensor T+ + of (tqu) along 
the shock wave. At a + = Oq", d+(f) (as well as <9+p) is discontinuous. T+ + will therefore 
have a delta function at <7q~. One easily finds 

T2 + (a+,a-) = MS(a+ - ff +)( ^J- } - 1) (eas) 

As the shock wave enters from J^((T _ = — oo), T+ + vanishes. As it moves in, P increases, 
and negative energy quantum fluctuations begin to accumulate along the shock wave. This 
energy approaches a constant up on which obeys 

T2 + +T f ++ = (tqf) 

at (oq,oo). Thus the / shock wave is gravitationally cloaked by a cloud of quantum 
fluctuations and arrives at as a "zero-energy bound state."* 



* This greatly strengthens the analogy made in [5,12] to the Schwinger model with a position- 
dependent mass. 
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This is in line with the conjecture made in [1] that the incoming state from evap- 
orates before forming a black hole, and arrives as a zero energy bound state on . It is 
tempting to speculate that this conjecture, though disproved for N > 24 [3,4], might be 
applicable to iV < 24. However it is not clear if this is consistent with the behavior of the 
apparent horizons. 

We hope to answer these questions in the near future. For the moment some further 
insight can be gained by investigating static black hole solutions of the type found for 
iV>24 in [6,7]. Following [7] this is best accomplished in terms of the variable 

s = A-V^-* 7 ^ (surb) 

which vanishes on the horizon, is spacelike and positive outside and timelike and negative 
inside. The equations of motion (pheq)- (rheq) for fields depending only on s become 

8P(s4> + <j>) = P'(\ 2 e 2p - 4s0 2 ), (speq) 

2P(sp + p) = -e- A(f> (X 2 e 2p - 4s0 2 ) , (sreq) 

while the constraint is 

e- 2 ^(4# - 20) = ^ (p 2 -p)-2((p - 0) 2 - p + $) + A {sceq ) 

where " " denotes differentiation with respect to s and t is a constant. Finiteness of 4> and 
p at the horizon gives constraints for initial data at s = 0: 

A 2 e 2p(O)-40(O) 



P(0) = 
0(0) = 



2P(0) 

A 2 e 2p(0)p/ (0) 



8P(0) 

i = 0. (hred) 

Since p(0) can be set to zero by a global coordinate transformation, there is a one parameter 
family of inequivalent solutions labelled by 0(0), or equivalently, the black hole mass. 
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The behavior of these solutions depends on whether 0(0) is less than or greater than the 
critical value, C = \ tn where P' changes sign. For large negative 0(0) (corresponding 
to large black holes) , the solutions will differ little outside the horizon from those found in 
[6,7]. Asymptotically the solution approaches the linear dilaton vacuum, but with a linearly 
divergent ADM mass corresponding to the infinite radiation bath required to stabilize the 
black hole against Hawking decay. As 0(0) approaches C the solutions begin to differ. 
This can be seen from the fact that at 0(0) = C the solution is exactly given by (in the 
gauge p(0) = 0) 

1 24 

^ = 2 in ~N " 

p= -£n(l + -). (dstr) 
a 

where a = 2A~ 2 (^ — 1). This corresponds to deSitter space filled with Hawking radiation. 
We presume that as 0(0) approaches C , there is a growing deSitter like region outside 
the black hole. At 0(0) = C , this region engulfs the entire spacetime and the black hole 
horizon becomes a deSitter horizon. 

One expects that a large slowly evaporating black hole is approximated within some 
region by these static solutions with a slowly increasing 0(0). However one should not 
conclude from the above that the endpoint of an apparent black hole formed by a massive 
incoming shock wave is deSitter space. In that situation, the spacetime is always asymp- 
totic to the linear dilaton vacuum with a finite mass. One possibility is that it looks like 
deSitter space within some region which then decays back to the linear dilaton vacuum. 

Inside the horizon, the ghost modifications have a crucial effect even for 0(0) << C . 
When N > 24, it was shown that increases (now in a timelike direction) until a zero of 
P and a singularity is reached. The resulting spacetime has a causal structure identical to 
that of the static classical black hole solution. However for N < 24 there are no zeros of 
P. In this case one finds from (speq) that for < C both and s0 are initially increasing 
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as before. will then inevitably cross C , at which point s0 must start to decrease, and 
"bounces" off of some maximum value rather than becoming singular. This behavior can 
be understood from a linearized analysis for small 0(0). To leading order the p equation 
is independent of and yields the deSitter solution (dstr). The linearization of (pheq) is 
then 

$<P + <t> = - f ^ ,o 0- (phln) 
(s + a) z 

Define the new timelike variable 



u = £n[ v ^ + v " ] (udet) 

\/OL ~ 



which runs from zero at the horizon to plus infinity at future timelike infinity (s = —a). 
The linearized equation then becomes 

d^4> + coth.ud u 4> + 20 = 0, (uphg) 

while the boundary condition (hred) implies d u (f> vanishes at the horizon. As the coefficient 
of the first derivative term is positive for u > 0, this is a damped harmonic oscillator. Thus 
excursions of (f> are damped inside the horizon and the linearized approximation does not 
break down. We therefore conclude that, at least for small 0(0), the interior of the black 
hole is non-singular and asymptotic to deSitter space, as illustrated in Figure 2. 

Evidently this system is very resistant to singularity formation: even if small black 
hole is forced into existence by continuously pumping in energy from infinity, there is no 
singularity in its interior. 

For large 0(0) the equations are harder to analyze, but we expect similar behavior. 
Numerical work on this question is in progress [8]. 

Clearly this set of equations exhibits complex and unusual behavior that we do not yet 
fully understand, and which merits further investigation. Our preliminary investigations 
have failed to uncover any black hole type singularities, but their existence is certainly not 
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ruled out. We also do not know if the equations give a qualitatively correct description 
of TV < 24 matter fields coupled to dilaton gravity, because higher loop corrections could 
be important. What has been established, however, is that the nature of the black hole 
formation/evaporation process, including the singularity structure, depends qualitatively 
on the properties of the matter sector. It is an urgent problem to characterize the possible 
behaviors. 

It is intriguing that the outcome of two-dimensional gravitational collapse depends 
qualitatively on the matter content of the universe. Perhaps this will also turn out to be 
true in four dimensions, and lead to constraints on the spectrum of elementary particles. 
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